Short-Range Correlations in "'He Liquid and Small "'He Droplets 
described by the Unitary Correlation Operator Method 
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The Unitary Correlation Operator Method (UCOM) is employed to treat short-range corre- 
lations in both, homogeneous liquid and small droplets of bosonic ''He atoms. The dominating 
short-range correlations in these systems are described by an unitary transformation in the 
two-body relative coordinate, applied either to the many-body state or to the Hamiltonian and 
other operators. It is shown that the two-body correlated interaction can describe the binding 
energy of clusters of up to 6 atoms very well, the numerical effort consisting only in calculating 
one two-body matrix element with Gaussian single-particle states. The increasing density of 
bigger droplets requires the inclusion of correlation effects beyond the two-body order, which 
are successfully implemented by a density-dependent two-body correlator. With only one ad- 
justed parameter the binding energies and radii of larger droplets and the equation of state 
of the homogeneous *He liquid can be described quantitatively in a physically intuitive and 
numerically simple way. 

PACS numbers: 67.40.Db 61.20.Gy 36.40.-c 
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I. INTRODUCTION 

The quantum-mechanical description of ground state 
properties of interacting many-body systems is a long- 
standing challenge in theoretical physics. In general the 
problem is two-fold: The first task is to specify the in- 
teraction between the constituents in terms of a two- 
body potential. The second is to solve the quantum- 
mechanical many-body problem with that interaction. 

For many-body systems composed of Helium atoms the 
first task, the construction of an appropriate interatomic 
potential, is solved rather precisely. The HFDHE2 poten- 
tial of Aziz et al. W is able to reproduce experimental 
transport coefficients in a wide temperature range. A 
simpler and less precise description of the interaction is 
given by the Lennard- Jones potential [g[ || 



vir) = 4e [{a/ry^ - {a/r)"] 



(1) 



with e = 10.22K and a = 2.556A. Since the emphasis 
of the investigations presented in this paper is more on 
novel concepts, which remarkably simplify the treatment 
of the many-body problem, we will restrict ourselves to 
this simple and widely used parameterization. 

A common property of many realistic microscopic in- 
teractions — like van der Waals type potentials be- 
tween neutral atoms or the nucleon-nucleon interaction in 
atomic nuclei — is a strong short-range repulsion, which 
induces strong short-range correlations in the quantum 



many-body state. The repulsive core prevents the parti- 
cles to approach each other closer than the radius of this 
core such that the two-body density matrix is strongly 
depleted for particle distances smaller than the core ra- 
dius. This effect cannot be described in a mean-field like 
picture or in terms of a superposition of a finite number 
of symmetrized product states. 

In this work we treat these correlations with the Uni- 
tary Correlation Operator Method (UCOM), that has al- 
ready been applied successfully in the framework of the 
nuclear many-body problem [y. 

In the following section we introduce the general for- 
malism of the Unitary Correlation Operator Method and 
derive analytic expressions for correlated wave functions 
and the correlated Hamiltonian. In the third section we 
determine the optimal correlation operator for liquid *He 
from the two-body system and discuss the structure of 
the effective interaction generated by the unitary trans- 
formation. In section IV we investigate the homogeneous 
■^He liquid with that correlator and introduce a density- 
dependent two-body correlator that describes very well 
the influence of three- and more-body correlations. Fi- 
nally, we successfully apply the density-dependent corre- 
lator to describe the ground state of small '*Hc droplets. 
Considering the great numerical simplification — only a 
few one-dimensional integrals have to be calculated — 
compared to Monte-Carlo-type calculations ||, ||], the 
agreement of the results is excellent. 
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II. UNITARY CORRELATION OPERATOR 
METHOD (UCOM) 

A. Concept 

The goal of the Unitary Correlation Operator Method 
(UCOM) is to describe short-range interaction-induced 
correlations by a unitary transformation of an uncorre- 
lated many-body state with TV particles. The unitary 
operator C associated with this transformation can be 
written as an exponential of a hermitian generator G 



C = exp[-iG]. 



(2) 



The generator is chosen to be an irreducible two-body 
operator 



N 



G=j:a, 



(3) 



1<J 



which describes genuine two-body correlations only. In 
principle we could add a three-body part to account for 
genuine three-body correlations. For the moment we 
omit this step to develop the basic formalism on a simple 
level. 

The particular structure of the unitary transformation 
— which should reflect the properties of the short-range 
correlations — is described by the two-body generator g. 
The repulsive core of the interaction leads to a suppres- 
sion of the probability density within the range Tcoic of 
the core. Uncorrelated many-body states that are sym- 
metrized products of single-particle states lead, however, 
to a two-body density that does not show any depletion 
for relative distances r < Tcore between the pair of par- 
ticles. Therefore the unitary transformation should shift 
those particles that are in the forbidden region of the re- 
pulsion to larger distances. This is done by means of a 
distance-dependent radial shift in the relative coordinate 
of two-particles. The hermitian generator of this shift in 
two-body space reads 

g^-[s{r){f-q) + {q-f)s{r)], (4) 

where q = ^(pi — P2) is the operator of the relative 
momentum, f — X1—X2 is the relative coordinate, and r 
is the associated unit vector. The function s{r) describes 
the size of the shift as function of the distance of the 
particles. 



B. Correlated Wave Functions 

First we want to apply the unitary correlation opera- 
tor C to an uncorrelated many-body state 1^) out of a 
low-momentum model space. The correlated many-body 
state \^) is defined by 



1^) =C|5'). 



(5) 



Here and in the following all correlated quantities are 
marked by a tilde. 

In order to demonstrate the generic effect of the corre- 
lation operator we discuss first a two-body system with 
an uncorrelated two-body state lip). Using the definition 
of the correlation operator C and the generator g (Q) we 
can evaluate the correlated two-body wave function in 
coordinate representation explicitly 



{X,r\C IV-) =TZ-{r) {X,R-{r)f\i:). 



(6) 



The respective transformation with the hermitian adjoint 
correlation operator C^ leads to 



{X,r\C^\i')=n+ir) {X,R+{r)f\^p). 



(7) 



The correlation operator acts in the relative coordinate 
r of the two-body system only, the center of mass coor- 
dinate X = ^{xi + X2) is invariant. According to the 
construction of the the generator (Eh the unitary trans- 
formation corresponds to a radial shift in the relative 
two-body coordinate 



R-{r) f , r 



R+{r)r. 



(8) 



The transformed relative distances are given by the corre- 
lation functions i?_ (r) and R^ (r) , respectively. Due the 
unitarity of the correlation operator CC^ = 1 the cor- 
relation functions i?_ (r) and _R+ (r) are inverse to each 
other 



R±[RTir)]^r. 



(9) 



The conservation of the norm of the correlated two-body 
wave function is guaranteed by the metric factor TZ±{r) 
in (|) and (§ 



n±{r) 



R±{r) 



v/^iM, 



(10) 



where R±{r) denotes the derivative of R±{r). It cor- 
responds to the square root of the Jacobi determinant 
associated with the transformation of the normalization 
integral. 

The correlation functions R±{r) that describe the co- 
ordinate transformation are connected to the function 
s{r), which enters in the generator (|j), by the integral 
equation 



±1 = 



For small shift distances the approximate relation 

R±{r) w r ± s{r) 



(11) 



(12) 



holds. 

Since all correlated quantities can be expressed directly 
in terms of the correlation functions R^{r) and R+{r) 
we will use one of these as basic quantity to describe the 
detailed structure of the short-range correlations. The 
shift function s(r) is only used for formal considerations 
and never specified explicitly. 



C. Correlated Operators 

Complementary to the definition of a correlated state 
(0) we can use the correlation operator to define corre- 
lated operators. For an arbitrary observable O the corre- 
lated operator O is given by the similarity transformation 



O^C^'OC. 



(13) 



It is evident that the formulations in terms of correlated 
states and correlated operators are equivalent when ex- 
pectation values or matrix elements are calculated 

(I'l O I*') = (*| C^OC I*') = (l-l O I*'). (14) 

Thus we can choose the formulation that is technically or 
intuitively better suited for the specific application under 
consideration. In most cases the formulation in terms of 
correlated operators is easier to apply in the many-body 
system. Moreover it provides a systemati c app roximation 



II D. 



scheme that will be discussed in section 

As for the correlated wave function we first discuss 
specific correlated operators for a two-body system. The 
generalization to many-body systems will be done in the 
following sections. 



1 Correlated Two-Body Potential 

The simplest operator of interest is the local two-body 
potential v{r) that depends on the relative distance r 
only. To evaluate the correlated potential we use the 
results on the correlated two-body wave function to de- 
termine a general two-body matrix element of the form 



(^|Ctw(r)C|V/) 

(d^r d^X (i/.| C^ \X, f) v{r) (X , r\ C [V-') 

U^rd^X (V^l^.f) v[R+{r)] (X,f|V/) 
{^\v[R+{r)]\^'). 



(15) 



In going from the second to the third line we use equa- 
tions (|6|) and (^ and substitute the integration variable r 
by i?_ (r) f. Since the two-body states are arbitrary this 
relation is valid on the operator level 



w(r) = Ctz;(r)C = w[i?+(r) 



(16) 



Thus the correlated two-body potential v{r) is given by 
the original potential with transformed radial coordinate 

r H^ R^{r). 



2 Correlated Kinetic Energy 

The correlated operator of the kinetic energy in the 
two-body system is slightly more complicated due to its 



momentum dependence. First we decompose the opera- 
tor of the two-body kinetic energy T2 in a center of mass 
contribution tcm and an relative part tioi 



*iol 



2M 



-P' 



2fi 



(17) 



The operator of the center of mass momentum P = 
Pi +P2 commutes with the correlation operator, i.e., the 
operator of the center of mass kinetic energy is invari- 
ant under the unitary transformation. The correlation 
operator acts on the relative part t^^i of the kinetic en- 
ergy only. A similar calculation as for the local potential 
leads to the following structure of the correlated relative 
kinetic energy B 



CH.ci C = trci + in + ir + u(r) 



(18) 



Besides tj-d an additional momentum dependent term ap- 
pears in the correlated relative kinetic energy, which can 
be formulated in terms of a tensorial effective mass cor- 
rection. It is conveniently split into a correction for the 
radial part tr and a different correction for the angular 
component t^ 



tn 



1 



l^ 



2jln{r) r2 



ir-i^-r)^f-cr). (19) 



where I ~ f x q is the operator of the relative two-body 
angular momentum. The distance dependent effective 
mass corrections fj-nif) and fj-rif) depend on the corre- 
lation function R+{r) only 



A* 



- 1 



M 



Mr) Rlir) ^' Mr) R?ir) 



- 1. 



(20) 



In addition to these momentum dependent terms a local 
contribution u{r) appears in ( p^ 



u{r) 



1 



1 



2fiR 



'21 



. KJr) 

'rR'M) 



R'Ur)\^ 1 R'^ir) 



R'M)J 2 R'M) 



(21) 



This additional two-body potential is also determined by 
the correlation function R+(r) and its derivatives. 



D. Cluster Expansion and Two-Body 
Approximation 

In the preceding section we evaluated correlated oper- 
ators explicitly in the two-body system. However, corre- 
lated operators in the many-body system contain addi- 
tional terms. 

The unitary transformation of an operator, e.g., an 
one-body operator like the kinetic energy or a two-body 



operator like the potential, generates a correlated opera- 
tor, which contains irreducible contributions to all parti- 
cle numbers 



o = c^oc 



(22) 



We use the notation O''"! for the irreducible fc-body part 
of the correlated operator O. This decomposition ac- 
cording to irreducible particle number is called cluster 
expansion y, 0]. 

The cluster expansion gains physical meaning by the 
so called cluster decomposition principle H]: Two local- 
ized subsystems, which are separated beyond the range 
of interactions, are independent of each other. There- 
fore the state of the total system decouples into a direct 
product of the states of the two subsystems. This implies 
an analogous decomposition property of the correlation 
operator and of correlated operators. 

Because of the cluster decomposition principle the clus- 
ter expansion is a natural starting point for approxima- 
tions of the correlated operators (p2). For a selected 
particle there will be a certain number of other parti- 
cles within the range of the correlation depending on the 
density. The number of particles in this cluster gives 
the maximum order of the cluster expansion that con- 
tributes. We can truncate the cluster expansion at low 
orders if the range of the correlations is sufficiently small 
compared to the average distance of the particles. 

The simplest nontrivial approximation results from the 
truncation of the cluster expansion beyond two-body or- 
der. The correlated operator in two-body approximation 
reads 



6^''' = 6^^^ + o^'-l 



(23) 



This approximation requires that the system is suffi- 
ciently dilute such that contributions of the three-body 
order of the cluster expansion are small. Or in other 
words, the probability for three and more particles to be 
in the range of the repulsive core simultaneously has to 
be small. The technical advantage of the two-body ap- 
proximation is that closed analytic expressions for the 
correlated operators can be deduced by just consider- 
ing the two-body system. The correlated Hamiltonian in 
two-body approximation is of the form 



H 



C2 



f[l] 
N 



N 



2] 



N 



Ef'+Ef+E*. 



r.[2] 



(24) 



i<3 



i<] 



The one-body part of the correlated kinetic energy is sim- 
ply the uncorrelated kinetic energy operator. The two- 
body part consists of the effective mass corrections ( [T9| ) 
and the additional local potential (^l]) as determined in 
the previous section 



f[2l =tj2+t;. +7-i(r). 



(25) 



The correlated potential has only a two-body contribu- 
tion that is given by (|l6| ) 



,[2] 



v(r). 



(26) 



In order to get a rough measure for the validity of the 
two-body approximation we define a smallness parameter 



K = pVc 



(27) 



as a product of the density p of the system and a typical 
volume in which the correlations between a pair of parti- 
cles change their relative wave function. The correlation 
volume Yc is defined by the norm of the defect wave func- 
tion, i.e., the difference between the uncorrelated uniform 
wave function (r (/)o) = 1 and its correlated companion 



Vc 



/"dV [7^_(r) -\f = fd^r [n+{r) - if 



(28) 



The smallness parameter k is a measure for the probabil- 
ity to find a third particle within the volume where the 
correlations between two particles change their relative 
wave function significantly. The two-body approxima- 
tion is valid only if this probability is small such that 
three-body correlations are negligible. We have shown 
for different physical systems [Q, || that the relative con- 
tribution of the three-body order to the energy exceeds 
10% if the smallness parameter reaches a value of typi- 
cally K = 0.3. 



E. Many-Body Correlations 

If the two-body approximation is not sufficient one 
can include higher orders of the cluster expansion succes- 
sively. However, already the calculation of a three-body 
correlated wave function starting from the many-body 
correlation operator 



II B is not 



in analogy to section 
practicable. Therefore we reverse the procedure and start 
from a general many-body coordinate transformation — 
in analogy to the correlation function R±{r) — and de- 
termine correlated many-body wave functions and opera- 
tors. In a second step we connect the explicit structure of 
the coordinate transformation with the many-body cor- 
relation operator, at least in an approximate way. 



1 Many-Body Coordinate Transformation 

Consider a iV-body system with a collective coordinate 
vector X = (xi, a?2, . . . ,a;jv)- The short-range correla- 
tions are described by a A'^-body coordinate transforma- 
tion 



X ^ X' = X^{X). 



(29) 



The transformation function X-{X) corresponds to the 
correlation function i?_ (r) in the two-body case. Similar 



to R+ir) we define an inverse transformation function 
X+{X) by 



X+[X^{X)] =x. 



(30) UiX) 



The correlated A^-body wave function (Xl^) is defined 
via the coordinate transformation (P£ 



(31) 



{X|*) ^V^{X) {X^{X)\^) 



with a metric factor T?-(X) that ensures norm conserva- 
tion. The metric factor is given by the square root of the 
Jacobi determinant of the transformation 



P4X) = ||J_(X)||V^ 



with the matrix elements 



[J-iX)] 



dX[ d[X_{X)], 
dXi 



dX, 



(32) 



(33) 



Formally we can construct a unitary many-body correla- 
tion operator C that generates this particular coordinate 
transformation 



/' 



C= /d^^X \X) V_{X) {X_{X)\ 



(34) 



The check of the unitarity relation C^C = 1 with this 
definition is straightforward. 

Using this formulation wc can evaluate many-body cor- 
related operators in the same way as shown in section |IIC| 
for the two-body approximation. We will only present 
selected results that are needed for the following investi- 
gations. The full formalism is discussed in m. 

A general local potential V{X), which depends on the 
coordinates X = {xi,X2, ■ ■ ■ ,xn) of all A^ particles, is 
transformed as in the two-body case 



CWiX)C = V[X+iX)l 



(35) 



i.e., the correlated potential is given by the original po- 
tential with transformed coordinate dependence. 

The expression for the TV-body correlated kinetic en- 
ergy operator C^TC is more involved 



3W 



C^TC=J2 Pj W{X)],k Pk + U{X) 



(36) 



j,fc=i 



where P = (pi,p2, • • • ,pW) denotes the collective mo- 
mentum vectors of all N particles. As in the two-body 
case ( p^ the correlated kinetic energy contains an effec- 
tive mass tensor and an additional local potential. The 
effective mass tensor is given by the square of the Jacobi 
matrix 



[M{X)],k 



1 
2m 



3N 



Y,[J-{x)h[J-{x)], 



(37) 



The additional local potential U{X) given by 

& 



1 
2m 



V-{X 



Wr\ Y.7^^'^-^^') 



i=l 



dX. 



X'=X+iX) 

(38) 



With these expressions it is in principle possible to calcu- 
late expectation values of the correlated Hamiltonian up 
to arbitrary order in the cluster expansion. In practice 
all extensions beyond two-body approximation are very 
costly, thus only the three-body order will be considered. 



2 Three- 



Correlations 



We will use the general formulation of the many-body 
coordinate transformation to estimate the contribution of 
the three-body order of the cluster expansion. Since we 
cannot derive the many-body coordinate transformation 
directly from the correlation operator we will construct 
the transformation by generalization of the known one in 
the two-body system. 

It is useful to reformulate the two-body coordinate 
transformation discussed in section |IIB| in terms of the 
general transformation function X\_ ' {X) , where the up- 
per index indicates the number of particles involved. The 
transformation (g|) of the relative coordinate with the 
correlation function R±{r) is equivalent to the following 
transformation of the single-particle coordinates 



A'i 



Xi 



Xi +A±{fi2) 



x\ 



X2 



A±(f2i 



X2 



with a shift vector 



A±(f) = i[i?±(r)-r]f. 



(39) 



(40) 



This notation reveals the intuitive picture behind the de- 
scription of short-range correlations by means of a coordi- 
nate transformation. As a consequence of the short-range 
repulsion the particles arc displaced along their connect- 
ing axis by the shift vector A±(r). 

The two-body transformation (B9|) can be readily gen- 
eralized to describe three-body correlations. The corre- 
spending three-body coordinate transformation Xj^ {X) 
is given by 



X 



(3) 



Xi 



Xl 



+ A+(fi2) + A+(fi3) 



A-i 



X2 



A+(r2i) + 



X2 



+ A+(f23) 



(41) 



A-V 



X3 



A+(r3i) + A+(r32) 



3^3 



Thus the coordinate of the first particle is transformed 
with respect to the second and the third particle simulta- 
neously. As a necessary prerequisite this transformation 
obeys the cluster decomposition principle. If one of the 



three particles is separated beyond the range of the corre- 
lations, i.e., to distances where A-|-(r) vanishes, then the 
transformation reduces to the two-body transformation 
(p9) for the remaining pair. 

We will use this three-body transformation to calcu- 
late the local three-body contributions of the correlated 
Hamiltonian, i.e., the correlated potential and the local 
part of the correlated kinetic energy. The irreducible 
three-body part ■ul'^1 of the correlated potential is given 
by 



i}!^! ^ C"^ {vi2 + vi:i + V2z) C 



-(*[1 



~[2] 



^23 



), 



(42) 



~ [21 

where v.ij is the uncorrclatcd two-body potential and v\- 
is the two-body part of the correlated potential. The first 
term describes the fully correlated two-body potential 
in a three-body system, which is given by a three-body 
coordinate transformation according to (Bq) . The second 
term is the sum of the two-body correlated potentials 
(EQ) . An analogous expression results for the local three- 
body part of the correlated kinetic energy til^l, which 
involves (|3|) and <^. 



III. OPTIMAL CORRELATION FUNCTION 

The only input needed to evaluate correlated quanti- 
ties explicitly is the correlation function R±{r). Since 
R± (r) reflects the properties of the correlations induced 
by the two-body interaction it can be extracted from the 
interacting two-body system. 



A. Mapping of the E = Scattering Solution 

One method to construct an optimal correlation func- 
tion is based on the exact solution of the two-body 
Schrodinger equation. The exact relative wave function 
4>e-^{f) contains all information on interaction induced 
two-body correlations. For the case of energy E = Q 
and relative angular momentum I = the solid line in 
Figure ^ shows the radial wave function (solid line) for 
the Lennard- Jones potential (|l|) (thin dotted line). The 
short-range part is dominated by a huge correlation hole, 
i.e., a region where the strong repulsive core of the po- 
tential enforces a vanishing wave function. 

We can extract a correlation function that describes 
these short-range correlations by requiring that the ex- 
act solution is reproduced by a correlated ansatz wave 
function 

(r|C|</)o) = (f|0ex) or (f|0o) = (^|Ct|,/)e,) (43) 

for all r inside some maximum radius A. Thus the op- 
timal correlation operator maps the short-range part of 
the uncorrelated wave function 00 (^) = (m'^o) onto the 



exact solution 4)cxij) = (r 0cx)- In coordinate represen- 
tation this condition leads with (0) and ( [lO|) to 

rMr) = i?+(r)y^i?V(r-)</>cx(r), (44) 

which can be written as an implicit integral equation for 
the optimal correlation function R+{r) [15] 






MO 



1 2 



c[i?+(e)] 



(45) 



The uncorrelated wave function (/jq (r) should be chosen in 
accordance with the uncorrelated many-body trial state. 
For example, if a condensate of N uncorrelated bosons 
is described by a iV-fold product of a Gaussian-shaped 
single-particle wave functions the relative wave function 
for each pair of bosons is again of Gaussian shape. 

The choice of the uncorrelated wave function (po (r) de- 
termines, which structures of the exact solution are gen- 
erated by the correlation operator and which have to be 
described by the ansatz state. In order to isolate short- 
range correlations the uncorrelated state should be able 
to describe the long-range behavior of the exact solution. 

According to this requirement we construct a hybrid 
wave function with a long-range behavior given by the 
exact zero-energy solution outside the range A of the in- 
teraction: 



(t>cxif) = 1 — a/r for r > A, 



(46) 



where a is the s-wave scattering length of the interac- 
tion. For two ''He atoms interacting via the Lennard- 
Jones potential (|l|) we get a = —171.3 A. The short-range 
part of the ansatz state is described by a Gaussian. At 
some radius % both functions are matched with continu- 
ous derivative. The final form of the uncorrelated wave 
function reads 



Mr) 



-^Mw^i^-^)] ■^r<x 



1 — a/r 



(47) 



■,r>x 



The matching radius x is chosen such that the correlation 
function R^{r)~r determined from (E^) vanishes at some 
finite radius A. This allows a natural separation of long- 
and short-raiigc correlations in the two-body system. 

In Figure U the hybrid ansatz with matching radius 
X = 7.83 A is shown in comparison with the exact E = 
solution. Both wave functions agree nicely at large dis- 
tances. At small distances the exact solution shows the 
pronounced depletion within the core and an enhance- 
ment in the attractive region of the potential. These 
short-range correlations are not contained in the uncor- 
related wave function (dashed curve). Instead the uni- 
tary correlation operator C, which maps the uncorrelated 
wave function onto the exact eigenfunction (solid curve) , 
takes over the task. As discussed in section [I C this has 



the great advantage to allow the definition of correlated 
interactions that act between simple uncorrelated many- 
body states like product states. 
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FIG. 1: Exact relative wave function <j)cx{r) for the E — 
scattering problem (solid line) in comparison with the hybrid 
ansatz (ftoir) for ^ = 7.83 A (dashed line). In addition the ra- 
dial dependence of the Lennard- Jones potential v{r) is shown 
(thin dotted line). 



FIG. 2: Optimal Correlation Function R+{r) — r resulting 
from the integral equation (Uq) (solid line) and the fitted pa- 
rameterization (thin solid line). For comparison the param- 
eterization with parameters resulting from energy minimiza- 
tion (thin dotted line) is depicted. 



The solution of the integral equation (|4^) for the cor- 
relation function R+{r) using these two wave functions 
is shown by the dashed line in Figure g. The quantity 
i?_i_ (r) — r corresponds to the net distance by which two 
particles with an initial separation r are shifted away 
from each other. The shape of the curve shows that 
particles with distances smaller than the core radius 
(''core ~ 2. 5 A) are shifted out of the core. If the initial 
distance is already larger than the core radius, then the 
shift decreases rapidly and finally vanishes at r « 6.6A. 
The small shift for larger radii (r > 6.6A) originates from 
the minimal deviations of exact wave function and hybrid 
ansatz in this range. However these are long range cor- 
relations, which should not be described by the unitary 
correlation operator. Therefore we will set the function 
R+{r) — r to zero for r > 6. 6 A in order to define the 
short-range correlation function. 

This intrinsic separation between long- and short-range 
correlations works only if the potential has a sufficiently 
strong attractive region in addition to the repulsive core. 
For this class of interactions the exact two-body wave 
function is depleted inside the core and enhanced in the 
attractive region. Since the unitary correlation operator 
conserves the norm of the wave function it can generate 
this kind of structure very easily: The amplitude shifted 
out of the core is placed in the attractive region. 

For purely repulsive potentials the situation is differ- 
ent. The exact wave function of a scattering state is de- 
pleted inside the core and the excess probability is spread 
over a large volume. Therefore the unitary mapping leads 
to a correlator of large range. In order to stay within the 
two-body approximation one has to restrict the range of 
the correlations and compensate by improving the uncor- 
related trial states. 



B. Energy Minimization 

An alternative method to determine the correlation 
function R+{r) emerges from the variational principle. 
If we think of the energy expectation value to be calcu- 
lated with correlated states and the uncorrelated Hamil- 
tonian then the correlation function introduces additional 
degrees of freedom into the trial state. By choosing a 
suitable parameterization for i?+ (r) with few variational 
parameters we can minimize the correlated energy expec- 
tation value and obtain an optimal correlation function. 

According to the structure of the correlation function 
determined by the mapping procedure described in the 
previous section we use the following parameterization 

R+{r)=r + a (r//?)'' exp[- exp(r//3)] (48) 

with three free parameters. The parameter (3 determines 
the range of the correlation function, a is related the 
maximum shift distance and thus the radius of the core, 
and r] influences the slope of R+{r) at small radii, i.e., the 
"hardness" of the core. We used this parameterization 
successfully for the description of the structure of atomic 
nuclei and nuclear matter 0, M. 

In order to describe a homogeneous ^Hc liquid we 
choose a constant wave function as two-body trial state 
(f|0o) = 1. In that case only the local potentials v{r) 
and u{r) contribute to the expectation value of the cor- 
related Hamiltonian 



£; = (0o| -ff |0o) = (00 1 i{r) + u{r) |(/)o) 
= \d?r [v{r) +u{r)]. 



(49) 
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FIG. 3: Local components of the correlated Hamiltonian: un- 
correlated Lennard- Jones potential r'^v{r) (thin dotted line), 
correlated potential r^v{r) (dashed line) and sum of corre- 
lated potential and local part of the correlated kinetic energy 
r^{v{r) + u{r)) (solid line). 



The expectation values of all momentum-dependent con- 
tributions t,:c\,ta,tr in (nm vanish. By inserting the pa- 
rameterization of the correlation function into (|l^) and 
( pl| ) we can express the integrand as function of the three 
variational parameters a, (3 and rj. The integration as 
well as the minimization can easily be done numerically. 
We obtain a unique set of parameters for the optimal 
correlation function 



6.267 A, /3 = 3.520A, 77 = 0.052. 



(50) 



The solid line in Figure y] shows the resulting correla- 
tion function. It is in very good agreement with the cor- 
relation function determined by the mapping procedure 
(dashed line). This demonstrates that the parameteriza- 
tion used in the variation is well suited for the correlation 
function. 

For numerical simplicity we will use the parameterized 
correlation function (Bq) with parameters according to 
(pQ) for the following investigations. 



C. Correlated and Effective Interaction 

With the explicit form of the correlation function we 
can study the structure of the correlated Hamiltonian in 
detail. As shown in section II C the unitary transforma- 



tion generates a momentum-dependent effective interac- 
tion. 

The local part of this interaction consists of the cor- 
related interaction potential v{r) (11^) and a local contri- 
bution u{r) (|2^) from the correlated kinetic energy op- 
erator. Figure M compares the radial dependencies of 
the uncorrelated Lennard-Jones potential, the correlated 



FIG. 4: Radial dependencies of the effective mass corrections: 
angular effective mass correction n/jln (r) (solid line) and ra- 
dial effective mass correction ^/jlr(r) (dashed line). 



Lennard-Jones potential and the sum of all local con- 
tributions. The unitary transformation of the potential 
shifts the repulsive core to smaller radii. This leads to a 
substantial reduction of the strength of the core, only a 
moderate repulsion at very short ranges remains. If we 
add the local part of the correlated kinetic energy the 
repulsive part vanishes completely; the local component 
of the correlated Hamiltonian is purely attractive. 

Actually, it can be shown analytically that mapping of 
an uncorrelated wave function, which is constant within 
the range of the repulsive core, onto the exact E = 
solution results in a local part v{r)+u{r) of the correlated 
Hamiltonian that vanishes identically inside the core [Q. 
The small contributions of the correlated local potentials 
within the core range that show up in Figure 
from the use of the restricted parameterization ( ^ 
the correlation function and the slight variation of the 
trial state. 

The momentum-dependent parts of the correlated 
Hamiltonian. which are formulated in terms of distance- 
dependent effective mass corrections (pO|), are shown in 
Figure |[ The radial effective mass correction (dashed 
line) generates a repulsion if the uncorrelated relative 
wave function has no n- vanishing radial derivatives. The 
angular effective mass correction (solid line) generates 
an attraction if the uncorrelated wave function has non- 
vanishing relative angular momentum. In general the 
contribution of these components is small compared to 
the total energy expectation value as the uncorrelated 
states have small gradients. 

Formally we can combine all the two-body terms of 
the correlated Hamiltonian to a momentum-dependent 
effective interaction 

w^'^'^ = v{r) + u{r) + in + ir (51) 

with a purely attractive local part given by (Hq) and (|2l|) 




and a momentum dependence given by (|lj). With this 
effective interaction the correlated Hamiltonian in two- 
body approximation has the standard form 



H 



C2 



N 

E 



~ 2 
W ■ 



(52) 



The next orders of the cluster expansion contribute ad- 
ditional three-body, four-body, and higher order effective 
interactions. 

The tamed effective two-body interaction i&I^l has a 
different operator structure and completely different ra- 
dial dependencies than the uncorrelated potential v{r). 
However, both interactions are phase-shift equivalent, 
i.e., in a two-body scattering process both interactions 
generate identical phase-shifts for any collision energy, 
by construction. This generic property of the effective 
interaction results from the finite range of the unitary 
transformation, which does not influence the asymptotic 
behavior of the scattering solutions. It is completely in- 
dependent of the shape of the correlation function at 
short ranges. In that way the unitary correlator provides 
a recipe to generate an infinite manifold of phase-shift 
equivalent interactions. This aspect is of special inter- 
est in cases where the interaction is determined only on 
the basis of two-body scattering data — like the nucleon- 
nucleon interaction. 



IV. HOMOGENEOUS *HE-LIQUID 

As a first application of the UCOM formalism we want 
to investigate the equation of state, i.e., the energy per 
particle as function of the density, for a homogeneous "^He 
liquid at temperature T ~ OK. 

In this case all effects that the interaction has on the 
many-body state are described by the unitary correlation 
operator. The uncorrelated many-body state pl'o) is a 
direct product of identical constant one-body states |'0o) 



|*o) = IV'o)®---® IV'o)- 



(53) 



For the calculation of expectation values we assume a 
finite volume V containing N particles. Accordingly the 
constant one-body states read 



{x\^o)^V-'/'. 



(54) 



The limit N,V ^ oo a.t constant density p = N/V is 
performed in the final step of the calculation. 



A. Equation of State in Tavo- and Three-Body 
Approximation 

First we calculate the energy per particle with the cor- 
related Hamiltonian in two-body approximation 



HC2^j,^ym^^[2]^fl2] 



fpl. 



The energy per particle in two-body approximation is 
defined by the expectation value of this Hamiltonian cal- 
culated with the uncorrelated many-body state (m) 



£^^2 



ip) 



1 

N 



(*o|fl-^'|*o) 



(56) 



Since the uncorrelated state is a product of constant 
one-body states only the local components V and U 
contribute; the expectation values of all momentum- 
dependent terms vanish. A simple calculation for the 
expectation value of the local components leads to the 
following expression for the energy per particle in two- 
body approximation 



£'^^(p) = (C[^l+C[^l)p 
with constant coefficients given by 



(57) 



CF - I jd\ vir), Cf^ = 1 Jd'r u{r). (58) 

The integrals over the correlated potential v{r) given by 
( p^ ) and the local part of the correlated kinetic energy 
vjr) according to (Uw are evaluated numerically 



Cl,2l = -856.5 K A3 , Cpl = 40 1 . 5 K A^ 



(59) 



Obviously, for homogeneous Bose fluids at zero temper- 
ature the two-body approximation is not able to describe 
saturation, i.e., a minimum of the energy at some finite 
density. The energy in two-body approximation ( pT\ ) is 
always proportional to density and due to the negative 
sum of the coefficients (p9) drops with increasing density. 

The dashed line in Figure ^ shows this behavior in com- 
parison with results of a Green's Function Monte Carlo 
(GFMC) calculation Q, which are exact up to the sta- 
tistical errors of the Monte Carlo sampling. Despite the 
fact that the two-body approximation fails to generate 
a minimum, the predicted orders of magnitude for the 
energy are correct. This is remarkable, considering that 
the uncorrelated Lennard- Jones potential would give an 
infinite positive energy expectation value for all densities. 
The special unitary transformation that we apply tames 
the strong core of the potential and leads to a bound 
many-body system. 

We note that the two-body approximation does gener- 
ate saturation in other systems like homogeneous Fermi 
liquids. In these cases the mass corrections contribute, 
because of the Fermi motion, and generate a positive con- 
tribution to the energy that grows with a higher power 
in density than the attractive local terms. 

A straightforward step to go beyond the two-body ap- 
proximation is the inclusion of the three-body terms of 
the cluster expansion. In addition to the terms in (|55| ) 
the correlated Hamiltonian in three-body approximation 
contains three-body contributions from the correlated 
potential and the correlated kinetic energy 



(55) 



H 



C3 



JjC2_^yl3] _^lfl3] + fP]^ 



(60) 
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FIG. 5: Energy per particle for a homogeneous ''He liquid 
as function of density calculated in two-body approximation 
e (p) (dashed line) and in three-body approximation e (p) 
(dotted line). The vertical line marks the density at which 
the smallness parameter k reaches the value 0.3. The circles 
show the results of an exact GFMC calculation [p|; statistical 
errors are smaller than the size of the symbols. 



Again only the local three-body terms V'I'^1 and C/'^l con- 
tribute to the energy expectation value with the uncor- 
related states (^). The resulting energy per particle in 
three-body approximation reads 



£^^3 



(p) = (c|;^i + c|fi)p+(c|;^i + c|fi)p^. (61) 



The coefRcients are defined by the six-dimensional inte- 
grals 

Cf^ = i yd^xid^xs i)[3l {.xuX2,x^ = 0) (62) 

[31 

and the corresponding expression for Ci . The inte- 
gration is done with a standard Monte Carlo algorithm 
("VEGAS" from @) and yields 



Cf^ =(-2512±50)kA^ 
Cf} = (29739 ± 1100) KA^ 



(63) 



The local three-body terms of the correlated Hamiltonian 
generate a positive contribution to the energy per par- 
ticle which is proportional to p^ . The dominant term is 
the local three-body contribution of the correlated kinetic 
energy. Together with the negative density-proportional 
contribution of the local two-body potentials the energy 
in three-body approximation ( |6l|) provides a minimum 
at finite density. 

The density dependence of the energy per particle in 
three-body approximation e'^^{p) is shown by the dotted 
curve in Figure ^ Obviously the three-body contribu- 
tion is not a small correction to the two-body approxima- 
tion for densities larger than 0.005 A"-^, which is only one 



fourth of the saturation density. This demonstrates that 
many-body correlations have a strong influence in these 
systems and that the two-body approximation alone is 
not sufficient to describe ground state properties. A sim- 
ilar conclusion can be drawn from the smallness param- 
eter K = p Vc ■ It reaches the phenomenological limit 
K = 0.3 for the validity of the two-body approximation 
at about 1/3 of the expected saturation density . 

However, the naive inclusion of the three-body terms 
of the cluster expansion of the correlated Hamiltonian 
does not improve the result beyond a density of about 
0.005 A^'^. Although we obtain saturation in three-body 
approximation, both, energy and density at the mini- 
mum of s'-^^{p) are substantially smaller than the re- 
sults of the GFMC calculation. This discrepancy has 
several reasons: Firstly, the optimal correlation function 
was determined in the two-body system, i.e., in two-body 
approximation. For a consistent treatment on the level 
of the three-body approximation the optimal correlation 
function should be determined by minimizing the energy 
in three-body approximation. Secondly, we expect that 
genuine three-body correlations play a very strong role 
|10| , [ll| ] . They are described by an irreducible three-body 
part YliiKjKkQijk in the generator (0) that was not in- 
cluded in the present treatment. Finally, in view of the 
large three-body contributions there is no good reason to 
neglect terms beyond three-body order. 

Each of the points mentioned can in principle be explic- 
itly included in the calculation. However, none of these 
corrections are expected to yield reasonably converged 
results in the particular case of the *He liquid. The den- 
sity range where the description should be applied is far 
above the range where the two-body approximation is 
applicable, i.e., where the smallness parameter k is be- 
low the typical value of 0.3. The successive inclusion of 
higher orders of the cluster expansion will extend this 
range but one would have to include very high orders — 
which arc beyond the numerical possibilities — to reach 
the expected saturation density. The only reasonable way 
to extend the treatment is a partial summation over all 
cluster orders. This is done in the framework of Jastrow 
correlations with the so called Hypernetted Chain (HNC) 
summation schemes [0, |l^, |ll|, ^|. However, this type of 
partial summation is not feasible for the unitary correla- 
tion operator because of its more complicated operator 
structure. 



B. Equation of State with Density-Dependent 
Correlation Functions 

We aim at an effective description of the higher or- 
ders of the cluster expansion that still allows a compact 
analytic formulation of correlated observables and does 
not require extensive numerical efforts. One general way 
to simulate the effect of higher cluster orders is to intro- 
duce density-dependent correlation functions on the level 
of the two-body approximation. 
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At low densities the correlator shifts a pair of particles 
in an optimal way from the repulsive into the attractive 
region. When the density grows other particles will be 
nearby and obstruct the correlation that was optimal for 
a free particle pair. To simulate this effect the range 
of the correlation function R+{r) should be more and 
more reduced for increasing density. This will reduce the 
binding as the particles will not be able any more to fully 
exploit the attractive region of the potential. 

In practice this is accomplished by scaling the parame- 
ters a and /9, which have the dimension of a length, in the 
parameterization (48) of the correlation function R+{r) 
with some factor £^{p) < 1 that depends on density. At 
very low densities we assume that the higher-order con- 
tributions are negligible such that £,{p — > 0) = 1. With 
growing density the effect of the higher cluster orders 
grows such that the scaling factor ^(p) has to drop. The 
most simple ansatz for the density dependence of ^(p) is 



C(p) = 1 - 7P, 



(64) 



with one free parameter 7. The density-dependent cor- 
relation function reads 



R+{r,p) = r + a^ip) 
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FIG. 6: Energy per particle for a homogeneous ^He liquid 
calculated in two-body approximation with a two-body op- 
timized correlation function i (p) (dashed line) and with 
density-dependent correlation function £p{p) (solid line). 
The vertical line marks the density at which the smallness 
parameter k reaches the value 0.3. The gray circles show the 
results of an exact GFMC calculation pi, the minimum is 
indicated by the black circle. 



Using this correlation function all components of the cor- 
related Hamiltonian — like correlated potential, local 
part of the correlated kinetic energy, and effective mass 
corrections — become explicitly density-dependent. The 
energy per particle is given by 



e^2 
^p 



(P) 



(Pr [vp{r,p) +Up(r,p)], 



(66) 



where the subscript p indicates that the density- 
dependent correlation function is used. 

Due to the explicit density-dependence of the corre- 
lated local potentials the equation of state is not propor- 
tional to p any more. From a Taylor expansion of (pq) 
around p = we obtain terms to all powers of the density. 
Thus all orders of the cluster expansion are represented 
by the density-dependent correlator in an effective way. 

The phenomcnological parameter 7 in the scaling fac- 
tor ( pil) is in general adjusted such that (£^) agrees with 
one experimental data point, or for some density with 
the result of a realistic calculation. We will use the re- 
sult of the GFMC calculation |] for the density p^^™^ 
and energy e^^'^^ at the saturation point and require 



xC2i CiYMC\ 
'-p iPsat / 



.GFMC 



(67) 



This ensures that the equation of state (|6g) calculated 
with the density-dependent correlator runs through that 
point. However, this does not imply that the minimum 
of (pq) _coincides with the GFMC minimum. Using the 



values |5| 
we obtain 



.GFMC 
Psat 



0.0219 A-3 and e 



3.696 A^. 



GFMC 



-6.848 K 



(68) 



We have chosen to fix the density dependence, i.e., 7, 
with a calculation instead of an experimental value, be- 
cause that calculation uses the same potential (|l|) and it 
is known that this simple Lennard- Jones potential does 
not reproduce exactly the experimental two-body data. 
Figure^ shows the energy per particle e?^(p) calcu- 
lated with the density-dependent correlation function. 
The phenomenological description of the higher order 
contributions by a density-dependent correlation func- 
tion indeed generates a minimum of the energy with 



Psat 



0.0212 A" 



?C2 

^p.sat 



-6.860 K. 



(69) 



The position of the minimum (indicated by the cross) as 
well as the shape of the equation of state around the min- 
imum agree very nicely the the GFMC result. As noted 
above this is not an immediate consequence of the adjust- 
ment of 7. The agreement shows that the phenomeno- 
logical density-dependence is well suited to describe the 
major aspects of many-body correlations on a numeri- 
cally very simple level. Actually, all calculations have 
been performed within a simple Mathematica notebook. 
Once the density-dependence is fixed the correlation 
function can be used to evaluate any correlated observ- 
able of interest. The analytic expressions for many corre- 
lated observables in two-body approximation are rather 
simple, but the evaluation of higher cluster orders be- 
comes complicated if not intractable. Thus the inclusion 
of the phenomenological treatment of higher order effects 
by means of the density-dependent correlator is a very 
useful tool. Interesting quantities are, e.g., the radial 
two-body distribution function g{r), the static structure 
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factor S{q) or the one-body momentum distribution n(k). 
The corresponding observables for fermionic Uquids were 
discussed elsewhere [||. 

In the foUowing section we will use the density- 
dependent correlation function to investigate the struc- 
ture of small droplets of ''He. This will demonstrate the 
universality of the concept. 



V. SMALL ^HE-DROPLETS 

In a next step we want to use the Unitary Correlation 
Operator Method to investigate the ground state prop- 
erties of small *He droplets. 

The many-body problem is treated in a simple varia- 
tional framework. The many-body trial state |^,-/V) is 
— like for the homogeneous liquid — assumed to be a di- 
rect product of identical one-body states. The one-body 
states are described by a Gaussian wave function 



{x\^) = (™)-3/4cxp[-fV(2a)] 



(70) 



with a common width parameter a as the only varia- 
tional degree of freedom. The Gaussian trial state is ap- 
plicable for small droplets with N < 70; larger droplets 
show a saturation of the central density, which cannot be 
modeled with this ansatz |1^, |lj| . An advantage of the 
Gaussian ansatz is that the two-body wave function 



{xi,X2\i^,^) = {X\^){r\^) 



(71) 



kinetic energy Tint, i-e., the total kinetic energy T re- 
duced by the center of mass contribution Tcm, can be 
expressed by the relative two-body kinetic energy irei 



N 



where 



tint 



N 



- T Tcm - j,r2^tr 



1 Y-. 



i ' 



i<j 



1 

2mN 



cl^ij 7 



N 



(74) 



(75) 



Due to the simple structure of the many-body trial state 
the expectation value of the correlated Hamiltonian (|7^) 
in the A^-body system can be expressed by that of the 
two-body system 



£;C2 ^ {^!^n\H>1^\^,N) 



N{N -I) 



(*,2|^£2|^^2). 



(76) 



Moreover we can use that the intrinsic Hamiltonian (|73| ) 
acts on the relative part {A (f?j of the two-body wave func- 
tion only. This leads to a simple form of the correlated 
energy expectation value of the A^-body system 



E' 



CI 



(Ar-l)(0|t,ci|0) 
A^(A^-l) 



(0|w(r)+u(r)+ir 10). 



(77) 



can be separated analytically into a center of mass com- 
ponent (A' $) and a relative wave function (n0) 



(r|0) 



(7ra/2) "' 'cxp 



-3/4 

(27ra)-V4 



cxp[ 



-r7(4a)]. 



(72) 



This simple trial state is able to describe mean-field 
like effects of the interaction, i.e., the spatial localization 
of the atoms. However, any correlation effect beyond 
the mean-field level has to be described by the unitary 
correlation operator. 



A. Binding Energy and rms-Radius in Two-Body 
Approximation 

First we use the two-body approximation with the cor- 
relation function (^) obtained from the two-body scat- 
tering solution to investigate the ground state properties 
of small droplets. The correlated intrinsic Hamiltonian 



structure 



jrG2 
-"int 



in two-body approximation has the 



H, 



C2 



tint 



y[2] _^ [/[: 



^[2] 



(73) 



where capital letters indicate summation over all parti- 



cles, like yP] = J2 



N ~[2] 
■ V ■ 
i<3 V ■ 



The uncorrelated intrinsic 



We used that the expectation value of the angular ef- 
fective mass correction io vanishes due to the spherical 
symmetry of the relative wave function (relative s-wave 
state). The two-body expectation value of the relative 
kinetic energy is given by 



(0|t,el|</.) = ^ /dV 



f) 2 



J__3_ 
2/i4a' 



(78) 



For the correlated two-body potential and the local part 
of the correlated kinetic energy we obtain 

((^|-D(r) +u(r) |0) = d^r [v{r) + u{r)] \{r\(l3)\'^ . (79) 

Finally, the expectation value of the radial effective mass 
correction reads 



(0| *r \4>) 



2fi J flr(r) 



(r|0) 



(80) 



For a fixed correlation function R+{r) the radial depen- 
dencies of the correlated potential v{r), the local part of 
the correlated kinetic energy u{r), and the radial effec- 
tive mass correction ^/ jlr{r) are known analytically (see 
The minimization of the energy E'-^^ as 



section 



[I CD 



function of the width parameter a as well as the calcu- 
lation of the integrals in equations (^) and (gO|) is done 
numerically. 
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Another interesting observable is the mean-square ra- 
dius of the droplet, defined by the expectation value of 
the operator 



Rr, 



1 ^ ^ - 2 1 



N 



N 

N^H'-'v (81) 



where Xcm = jfTlii ^i is the center of mass coordi- 
nate of the many-body system. The formulation in terms 
of the two-body operator r^ reveals a similar structure 
like for the intrinsic kinetic energy (|74| ) . The many-body 
product state under consideration allows a direct calcu- 
lation of the uncorrelated expectation value using the 
relative two-body wave function (|72[) 



= [(*,iV|H,„3|*,iV)] 



1/2 



N -1 

2N 



(0|r2|0) 



1/2 



N ■ 



2N 



■ 3a 



1 1/2 (82) 



The correlated ms-radius i?ms ™ two-body approxima- 
tion has a similar structure like the kinetic energy 



r: 



C2 



R„ 



mi 



(83) 



In addition to the uncorrelated mean-square radius Rms 
the unitary transformation generates a two-body contri- 
bution 



J. Cr, 



-^f:m-^^)-^.]- 



(84) 



i<j 



The expectation value of the correlated ms-radius (p2 
can again be expressed by the relative two-body wave 
function alone 



fC2 
ims 
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In order to study the quality of the two-body approxi- 
mation we first look at very small droplets with A^ < 10. 
For these droplet sizes an early variational Monte Carlo 
(VMC) calculation on the basis of the Lennard- Jones po- 
tential exists [pi . The authors use a Jastrow-type param- 
eterization of two-body correlations with a long-range 
form adjusted to the behavior of the exact two-body so- 
lution. The calculation of the many-body expectation 
value with these trial states involves a Monte Carlo inte- 
gration routine. The resulting energy expectation values 
show a statistical error of approximately 10%. 



Our results of the minimization of E^'^ (77) with re- 
spect to a using the system independent correlation func- 
tion given in equations ( |48| ) and (^) for A^ < 10 are 
summarized in Figure [7[ The upper panel shows the 
correlated energy per particle, the middle panel the cor- 
related rms-radius of the droplet and the lower panel the 
smallness parameter Koq versus the number of particles 
N. 
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FIG. 7: Energy per particle (upper panel), rms-Radius (mid- 
dle panel), and smallness parameter Kaq (lower panel) for 
small ''He droplets in two-body approximation with two-body 
optimized correlation function. The gray circles show results 
of a variational Monte Carlo calculation H , statistical errors 
for the energy are indicated by the gray bars. The lines should 
guide the eye. 



As for the homogeneous liquid the expectation values 
of the uncorrelated Lennard- Jones potential would di- 
verge. But the unitary transformation of the Hamilto- 
nian using the correlation function determined in the two- 
body system removes the strong short-range repulsion of 
the potential completely. With this tamed correlated po- 
tential and a simple product ansatz for the many-body 
wave function the energies of very loosely bound droplets 
actually agree very nicely for A^ < 6 with the VMC re- 
sult. For A^ > 6 the energy expectation value is below 
the VMC value and the difference grows with increas- 
ing particle number or density. Since expectation values 
are evaluated in two-body approximation, i.e., the full 
cluster expansion is truncated above two-body order, a 
basic property of the Ritz variational principle is lost: 
The variational minimum of the energy is not necessarily 
bounded from below by the exact energy eigenvalue. The 
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observed overbinding indicates that the higher orders of 
the cluster expansion have to give a sizeable positive con- 
tribution to raise the expectation value above the exact 
eigenvalue. 

In direct connection with the overbinding in two-body 
approximation the correlated rms-radii of energy mini- 
mized states are systematically smaller than the VMC 
result. Accordingly the densities are too high. In con- 
trast to the homogeneous liquid the two-body approxi- 
mation does not produce a pathological collapse towards 
high densities for small droplets. This is due to the pos- 
itive contributions of the kinetic energy and the effective 
mass correction that are absent in the homogeneous case. 

As a measure for the validity of the two-body approx- 
imation we used the smallness parameter k = pVc with 
correlation volume Vc defined by (Eq) . In order to spec- 
ify the smallness parameter for inhomogeneous systems 
we have to define a measure for the average density. 
One possible definition is the value of the density for 
a step-like density profile with same (uncorrelated) rms- 
radius as the original distribution. For identical Gaus- 
sian single-particle wave functions with width parameter 
a this equivalent density is given by 



3A^ 



3 1 



^ rms 



3/2 



An 



N 1 



5 iV- 1 a 



3/2 



(86) 



The smallness parameter n^q = PcqVc defined with this 
density is shown in the lower panel of Figure R as function 
of particle number. Like for the homogeneous liquid the 
energy in two-body approximation starts to deviate from 
the exact result if the smallness parameter Keq exceeds 
the value 0.3 (indicated by the dashed line). This again 
confirms that the two-body approximation is valid as long 
as the smallness parameter fulfills the condition n < 0.3. 



B. Energy and rms- Radius w^ith 
Density-Dependent Correlation Function 
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FIG. 8: Energy per particle (upper panel), rms- Radius 
(middle panel), and smallness parameter Kcq (lower panel) 
for small *He droplets in two-body approximation with the 
density-dependent correlation function fixed in the homoge- 
neous liquid. The gray circles show results of a variational 
Monte Carlo calculation Ipl. The lines should guide the eye. 



To account for the effect of many-body correlations 
in a simple but efficient way we employ the concept 
of dens ity-de pendent correlation functions introduced in 
section 



IV B for the homogeneous liquid. 



To apply the density-dependent correlator in an inho- 
mogeneous system we have to specify in which way the 
density-dependence should be evaluated. The most sim- 
ple approach is to insert the equivalent density of the 
droplet, as defined in (p6|), into the density-dependent 
correlation function. That means we neglect that the ef- 
fects from stronger many-body correlations in the center 
of the droplet do not exactly cancel the weaker ones at 
the surface. 

A more elaborate ansatz to account for the inhomo- 
geneity would be to implement the density-dependence 
in local density approximation, i.e., the correlation func- 
tion which acts on the relative coordinate of a particle 
pair depends on the local density at their center of mass 
position. This is not done here. 



The variational calculation of the ground state proper- 
ties of small "^He droplets including a correlation function 
that depends on the equivalent density (Pq ) is straight- 
forward. The parameterization ( pq ) of the correlation 
function i?_|-(r, pcq) is used with the parameters ( p0[) de- 
termined in the two-body system. The parameter of the 
density-dependent scaling function (|64| ) is taken from the 
investigations of the homogeneous liquid (68). The en- 
ergy expectation value in two-body approximation with 
density-dependent correlations is of the same form as for 
the density-independent correlation functions discussed 
in the preceding section. 

The results of the energy minimization with respect 
to the width parameter a are summarized in Figure |8[ 
The energy expectation values are in full agreement with 
the VMC calculation [gj for all particle numbers. The 
overbinding observed with a density-independent correla- 
tion function is completely compensated by the density- 
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TABLE I: Energy expectation values and rms-radii for small ^He droplets with different particle number A'' resulting from 
energy minimization in two-body approximation with two-body optimized and density-dependent correlation function. The 
width parameter a is given in units of A^, energies are given per particle in units of K, and the rms-radius in units of A. 





Two-Body Approximation 


Two 


-Body Approximation 


with Density-Dependent Correlator 


Ref. 


H 


N 


a 


£^7a 


f:C2 

' rms 


a 


Tint /A 


m/N 


Ul'^/N 


V^'^/N 


Ef^N 


fC2 

' rms,p 


E/N 


''"rms 


3 


32.33 


0.000 


5.71 


34.39 


0.176 


0.005 


0.238 


-0.418 


0.002 


5.88 


-0.027 


6.60 


4 


22.83 


-0.085 


5.12 


25.19 


0.271 


0.021 


0.551 


-0.919 


-0.076 


5.37 


-0.103 


5.75 


5 


19.81 


-0.199 


4.97 


22.70 


0.320 


0.038 


0.846 


-1.379 


-0.175 


5.29 


-0.200 


5.34 


6 


18.26 


-0.329 


4.90 


21.71 


0.349 


0.054 


1.122 


-1.807 


-0.282 


5.30 


-0.300 


5.17 


7 


17.31 


-0.468 


4.88 


21.32 


0.366 


0.067 


1.378 


-2.202 


-0.391 


5.34 


-0.407 


5.13 


8 


16.67 


-0.613 


4.87 


21.23 


0.375 


0.079 


1.614 


-2.567 


-0.500 


5.40 


-0.506 


5.20 


9 


16.20 


-0.761 


4.88 


21.31 


0.379 


0.088 


1.832 


-2.905 


-0.606 


5.46 


-0.589 


5.26 


10 


15.84 


-0.913 


4.89 


21.50 


0.380 


0.096 


2.033 


-3.218 


-0.709 


5.53 


-0.680 


5.32 


20 


14.42 


-2.500 


5.14 


25.13 


0.343 


0.118 


3.427 


-5.442 


-1.554 


6.20 






30 


14.01 


-4.125 


5.44 


29.15 


0.301 


0.109 


4.221 


-6.778 


-2.146 


6.75 






40 


13.81 


-5.760 


5.74 


32.97 


0.269 


0.098 


4.745 


-7.698 


-2.588 


7.20 






50 


13.69 


-7.399 


6.02 


36.59 


0.243 


0.085 


5.119 


-8.381 


-2.933 


7.60 






60 


13.61 


-9.040 


6.29 


40.02 


0.223 


0.076 


5.405 


-8.918 


-3.214 


7.95 






70 


13.56 


-10.683 


6.55 


43.29 


0.207 


0.068 


5.631 


-9.353 


-3.447 


8.26 







dependence. The correlated rms-radii of droplets with 
A^ > 6 arc consistently larger by 0.2 A than the VMC re- 
sult. Unfortunately the authors j^] give no estimate for 
the statistical or systematic errors for their procedure to 
obtain the rms-radii. For the smallest droplets our result 
for the rms-radii are below the VMC results. This may 
be caused by the restriction to Gaussian trial states ([72|), 
which are not able to describe the long-range exponential 
tail present in these extremely weakly bound systems. At 
the same time the insufficient trial state causes a slight 
underbinding of the smallest droplets. 

For completeness Figure g also shows the smallness pa- 
rameter Koq = PcqVc- Siucc both, thc correlation volume 
Vc and the density of the droplets are reduced by the 
density-dependent correlator the product stays small. In 
any case Keq has lost its meaning as a measure for three- 
body correlations as the density-dependence of the cor- 
relator already includes many-body correlations. 

Table | summarizes the results of the energy minimiza- 
tion with two-body optimized and density-dependent cor- 
relator for droplet sizes up to TV = 70. In addition to 
the correlated energy and the rms-radius the individual 
terms of the correlated Hamiltonian are shown. The at- 
tractive correlated potential and the repulsive local part 
of the correlated kinetic energy are the major contribu- 
tions and show a large cancellation. The intrinsic kinetic 
energy and the effective mass correction give rather small 
contributions. 



VI. SUMMARY AND CONCLUSIONS 

We presented a simple and elegant approach to de- 
scribe strong short-range correlations in a many-body 
system that result from the strong repulsive core of the 
two-body interaction. These correlations induce high- 
momentum components into the relative wave function 



that cannot be treated within a low-momentum model 
space. The correlation operator C generates thc short- 
range correlations by a unitary transformation in the rel- 
ative coordinate of each pair of particles. It shifts the 
particles out of the repulsive core of the two-body inter- 
action. 

The correlation operator is then used to define corre- 
lated observables, most notably the correlated Hamilto- 
nian. The unitary transformation yields a tamed effective 
two-body interaction, which consists of a purely attrac- 
tive local part and a weak momentum dependence. In 
addition three-body and higher order terms are gener- 
ated. For dilute systems, i.e., if the range of the correla- 
tions is small compared to the average distance between 
the particles, it is sufficient to include the terms up to 
two-body order. In high-density systems the higher or- 
der terms can be included explicitly or simulated by a 
density-dependent correlation function. 

In order to study the applicability of the method in 
many-body systems, where short-range correlations play 
a dominant role, we investigated the ground state proper- 
ties of homogeneous ^He liquids and small droplets. The 
correlation function is fixed by energy minimization in 
the two-body problem with a constant uncorrelated wave 
function. Alternatively the correlation function can be 
determined by mapping an uncorrelated wave function 
onto the exact energy E = two-body wave function. 

We use this correlation function to calculate the energy 
of the homogeneous "^He liquid as function of density. 
The uncorrelated many-body state is a direct product of 
constant one-body states, thus the influence of the inter- 
action on the state has to be described by the correlation 
operator alone. Already in two-body approximation the 
strong short-range repulsion of the Lennard-Jones poten- 
tial is tamed completely and we obtain a bound system. 
However, the two-body approximation does not generate 
saturation, i.e., the energy per particle decreases propor- 
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tional to density. To simulate the effect of many-body 
correlations we employ a density-dependent correlation 
function. The one parameter introduced by the density 
dependent scaling of the correlation function is fixed by 
adjusting the correlated energy for one value of the den- 
sity to the result of a GFMC calculation. The resulting 
equation of state shows a minimum and is in very good 
agreement with the predictions of the GFMC calculation 
over the whole density range. 

In a second step we use the density-independent as 
well as the density-dependent correlation function to per- 
form an ab initio calculation of the ground state structure 
of small ^He droplets in a simple variational framework. 
It turns out that the two-body approximation with the 
density-independent correlation function gives a good de- 
scription for small droplets with particle numbers iV < 6. 
For larger droplets an overbinding compared to a quasi 
exact variational Monte Carlo calculation appears. Due 
to the higher density in these droplets many-body corre- 
lations become relevant. If we use the density-dependent 
correlation function with one parameter adjusted in the 
homogeneous liquid then the energies and rms-radii are 
in very good agreement with the exact calculation for all 
particle numbers under consideration. 

We conclude that the Unitary Correlation Operator 
Method is a powerful and universal tool to handle short- 
range correlations in the many-body system. For suffi- 
ciently low densities one can use the two-body approxi- 
mation to construct the correlated Hamiltonian as well 
as other correlated observables in closed form. The re- 
sulting effective two-body interaction can be used for an 
ab initio description of the many-body system based on 
a low-momentum model space. A gross measure for the 
validity of the two-body approximation is the smallness 
parameter k (vTt), which should be below a phenomeno- 
logical limit of 0.3. 

For higher densities one can utilize density-dependent 
correlation functions to describe many-body correlations 
effectively within the two-body approximation. The ad- 
ditional parameter in the density-dependence can be 
fixed by one measured data point or by comparison with 
an exact calculation. Nevertheless, the resulting density- 
dependent correlator can then be used for a broad range 
of applications, e.g.. ab initio calculations of correlated 
observables like energies, densities or momentum distri- 
butions of different systems without further adjustments. 



The great advantage of the Unitary Correlation Operator 
Method with density-dependent correlation functions is 
that besides its transparent physics it requires only min- 
imal computational power compared to the huge com- 
putational resources needed for Monte-Carlo-type cal- 
culations. A set of Mathematica notebooks that illus- 
trate the application of the Unitary Correlation Oper- 
ator Method to ^He liquids can be found at the URL: 
http: //theory .gsi .de/~rroth/math/. 
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